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Reduction of chemical networks 
I. The case of molecular clouds 
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Abstract. We present a new method to analyse and reduce chemical networks and apply this technique to the 
chemistry in molecular clouds. Using the technique, we investigated the possibility of reducing the number of 
chemical reactions and species in the UMIST95 database simultaneously. In addition, we did the same reduction 
but with the "objective technique" in order to compare both methods. We found that it is possible to compute the 
abundance of carbon monoxide and fractional ionisation accurately with significantly reduced chemical networks 
in the case of pure gas-phase chemistry. For gas-grain chemistry involving surface reactions reduction is not worth- 
while. Compared to the "objective technique" our reduction method is more effective but more time-consuming 
as well. 
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1. Introduction 

Despite the significant growth of computer power over the 
last years, the modelling of the chemical evolution of pro- 
tostellar cores and protoplanetary discs is still a challeng- 
ing computational task. This is even more true if one aims 
at the simultaneous self-consistent modelling of the dy- 
namical and chemical evolution of these systems. 

The mathematical formulation of the chemical pro- 
cesses in space is realised by differential equations describ- 
ing chemical kinetics. They contain terms for all chemi- 
cal reactions expected to proceed efficiently in the inter- 
stellar medium (ISM). Along with abundances of all in- 
volved species, these terms depend on the rate coefficients 
that are summarised in special databases. The relevant 
astrochemical rate files, like the UMIST rate file or Ohio 
New Standard Database, comprise hundreds of species and 



thousands of reactions (Millar et al. 1991, 1997; Le Teuff ct 
al. |2000[ Terzieva k Herbst |1998[ Aikawa & Herbst |1999) ). 
Of course, one can only wonder if these databases contain 
all the information about reactions needed to accurately 
reproduce observed abundances of interstellar molecules. 
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The opposite question to ask is if all this information is 
always needed. 

In many astronomically interesting situations it is 
enough to follow the evolution of a limited number of 
species only. Intuitively, it seems to be apparent that to 
predict an abundance of N2H + , for example, one can af- 
ford computing Cio abundance with less accuracy or even 
ignore this molecule altogether. In most cases, though, 
modelling of the chemical evolution of a parcel of the in- 
terstellar gas under fixed conditions or in a limited prede- 
fined range of conditions is not very time-consuming, so 
one probably will not be eager to sacrifice any accuracy 
even for the sake of order of magnitude gain in compu- 
tational speed. However, in the case of dynamical mod- 
elling with changing density, UV intensity, temperature 
etc., a reduced chemical network may help to distinguish 
between feasible and non-feasible problems or at least to 
ease the computational burden if parallel processors are 
in use. One possible application for reduced chemical net- 
works would be the modelling of the evolution of magne- 
tised protostellar clouds or protoplanetary accretion discs, 
when it is necessary to compute the fractional ionisation 
self-consistently. Another example of dynamically impor- 
tant species is carbon monoxide (CO) which under certain 
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conditions can be the most important molecular coolant 
(see, for example, Glassgold & Langer 1973Q . 

Several authors have already used reduced chemical 
networks, based on the above arguments, however, with- 
out presenting accurate mathematical considerations (e.g., 
Gerola & Glassgold 1978; Howe, Hartquist, & Williams 
1994| ; El-Naw awy, Howe, & Millar |1997| ). Recently, Ruffle 
et al. (|2002D (hereafter RRPHH) and Rae et al. (|2002p 



Table 1. Model designations. 



(hereafter RBHPR) made a first attempt to reduce astro- 
chemical networks by so-called "objective reduction tech- 
niques", developed in combustion chemistry. They have 
shown that it is possible to quantify the reduction task by 
introducing a certain criterion that would allow in each 
particular case to indicate, abundances of which species 
must be followed in order to have a reasonable estimate for 
an abundance of a species under investigation. For static 
regions with physical conditions typical of diffuse clouds, 
they isolated the species set that contains tens instead of 
hundreds species. This reduction only introduces a factor 
of two uncertainty in predicted CO abundance. The same 
task has been solved for the fractional ionisation as well. 

We perform a similar analysis of the UMIST 95 chem- 
ical reaction database in order to find if it is possible 
to reduce the number of species and/or reactions in the 
chemical networks used to calculate ionisation degree and 
carbon monoxide abundance under conditions typical of 
molecular clouds. Unlike RRPHH and RBHPR, we partly 
focus on the medium that is more advanced toward the 
formation of a star. This situation is more complicated 
from the chemical point of view. First, pre-stellar clouds 
are often shielded from the interstellar UV radiation field. 
Dissociating photons, being an additional chemical factor, 
in effect simplify the chemical network, as they prevent 
many complex species from being formed. Reactions in 
grain surfaces are another factor that cannot be ignored 
in dense clouds. 

The paper is organised as follows. First, we describe 
our chemical model in Sect. 2. The species-based and 
reaction-based reduction methods are outlined in Sect. 3. 
Results of reductions obtained with both methods are pre- 
sented in Sect. 4 (ionisation degree) and Sect. 5 (CO). 
These results are discussed in Sect. 6. A conclusion fol- 
lows. 

2. Chemical model 

We take the rates for gas-p hase reactions from the 
UMIST 95 rate file (Millar et al. |1997[ ). In addition to gas- 
phase chemistry, gas-dust interaction and dust surface re- 
actions are considered in some of the models. The rates 
at which surface reactions proceed and the very applica- 
bility of the rate approximation to surface chemistry still 
remain a subject of debate. This is why many authors pre- 
fer to neglect surface reactions altogether (for instance, see 
Ai kawa fe Herbst |199S| ; Bcrgin & Langer |1997| Ber gin et 
al. 1995; Charnley et al. 2001). However, in the attempt 



Notation 


Meaning 


DIFF 


nu = 10 a cm~ a 


DENS 


n H = 10 7 cm~ 3 


GAS 


Gas-phase network (395 x 3864 a ) 


GAD 


GAS with accretion and desorption 




(543 x 4593) 


DUST 


GAD with surface reactions (543 x 4785) 


HM 


"High metals" 


LM 


"Low metals" 



^The network contains 395 species involved in 3864 reactions. 



actions seems to be unavoidable, as dust is now widely 
believed to play a crucial role in molecular cloud chem- 
istry, serving both as a sink for frozen-out species and as 
a catalyst for the s ynthesis of complex molecules (see, how- 
ever, Turner 2000). In this paper we seek not to reproduce 
any particular abundance pattern observed in molecular 
clouds, but rather to investigate the very possibility of re- 
duction when surface reactions are present. This is why 
we did not try to incorporate any new ideas about surface 
chemistry. 



2.1. Physical conditions and initial abundances 

We solve equations of chemical kinetics for the range of 
fixed physical conditions, listed in Table |l|, where model 
designations are given as well. Gas and dust tempera- 
tures are 10 K in all cases. No direct UV radiation is 
present, implying the conditions, typical of obscured re- 
gions, like cores of molecular clouds, where the visual ex- 
tinction Ay > 10 mag. Dust grains are assumed to be 
silicate-like spheres with a constant size ad = 0.1 /im, den- 
sity pd = 3 g cm" 3 , and 10 6 surface sites for adsorption 



(Hasegawa et al. 1992 ). Dust is supposed to constitute 1% 
of the gas density by mass. 

All abundances in the paper are given with respect 
to the number of hydrogen nuclei. We use "high metal" 
and "low metal" sets of initial abundances (e.g. Lee et 
al. |1998[ ) quoted in Table |. In the "high metal" set the 
standard solar elemental composition is corrected with a 
modest depletion of 2 for S and stronger depletions of 10 
for Na, 50 for Si, 60 for Mg, and 110 for Fe. The "low 
metal" values contain additional depletion factors of 100 
for each of these elements. The abundances of all elements 
but P and CI are taken from Aikawa et al. (1996). For 
P and CI we take values from Grevesse & Sauval (1998) 



to constrain the chemical network in a cold and dense 
molecular cloud core, taking into account gas-dust inter- 



and use the same depletion factors as for Fe. Hydrogen is 
assumed to be completely in molecular form initially. All 
other elements are present in atomic form at t = 0. 

Below details of the adopted model of gas-dust inter- 
action are provided. 
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Table 2. Adopted elemental abundances. 



Element 


"High metals" 


"Low metals" 


He 


9.75(-2) 


9.75(-2) 


C 


7.86(-5) 


7.86(-5) 


N 


2.47(-5) 


2.47(-5) 





1.80(-4) 


1.80(-4) 


s 


9.i4(-e) 


9.14(-8) 


Si 


9.74(-7) 


9.74(-9) 


Na 


2.25(-7) 


2.25(-9) 


Mg 


1.09(-6) 


1.09(-8) 


Fe 


2.74(-7) 


2.74(-9) 


P 


2.16(-8) 


2.16(-10) 


CI 


1.00(-7) 


1.00(-9) 



2.2. Gas-grain interactions 
2.2.1. Neutral accretion 

The accretion rate fc ac («) of the ith species is given by 
fcac(i) = Tralvth(i)n d S, 



where i>th(i) = ■\/8fcTg as /7rm i is the thermal velocity for 
the ith species, is the number of dust grains per unit 
volume, T gas is the gas temperature, is the atomic mass 
of the ith species, and k is the Boltzmann constant. The 
sticking probability S is assumed to be 0.3 for all neu- 
tral species except for H, He, and H2. The sticking co- 
efficient of atomic hydro gen is estimated from Eq. [3.7] 
by Hollenbach & McKee (1979). Sticking probabilities for 
helium and molecular hydrogen are assumed to be zero. 



2.2.2. Dissociative recombination on grain surfaces 

Following Umebayashi & Nakano ( 198C| ), we suppose that 



all grains have unit negative charge and thus can act as 
electron donors in the same reactions of dissociative re- 
combination that occur in the gas phase. The rate coeffi- 
cient for reactions of dissociative recombination on grain 
surfaces is given by 



fedr(i) = ^a 3 Tra d vth(i)ridCu 



where 



S 1 



1.671 ■ lo- 
aded 



(1) 



(2) 



with and Td being the d ust si ze and temperature, re- 
spectively (Rawlings et al. 1992| ). Here the summation 



is over different recombination channels; the quantity ay 
denotes the probability of a particular channel. Sticking 
probability S is 0.3 for all ions. Products of these reactions 
are assumed to return to the gas-phase immediately. 

2.2.3. Desorption from grain surfaces 

The two desorption processes taken into account in our 
model are thermal evaporation and cosmic ray induced 
desorption. 



The thermal evaporation rate kth(i) of the ith surface 
species is 



kth(i) = voexp(-T D (i)/T d ), 



(3) 



where Td is the dust temperature, kT^>(i) is the binding 
energy for physical adsorption of the ith species to the 
dust surface, vq is the characteristic vibrational frequency 
for the adsorbed species given (Eq. [3] from Hasegawa et 
al. gggg ). 

The cosmic ray desorption rate k CI (i) of the ith surface 
species can be estimated as 



k CI (i) = 3.16- 10~ 1! Voexp(-T D (i)/70K). 



(4) 



The values of Td for most species are taken from Hasegawa 
& Herbst ( |1993] ). There are some species in the UMIST95 
database for which desorption energies are not given in 
that paper. For these species, Td is interpolated from val- 
ues of chemically identical species or estimated as 



T D = 50 A 



(5) 



with Ai being the corresponding atomic number, which 
is a good approxim ation for species from Table 4 of 
Hasegawa & Herbst (|l993|) . 

2.2.4. Surface reactions 

The surface reactions as well as their rates are taken from 



Hasegawa et al. ( 1992 ) with the appropriate correction for 
different densities. Activation energy of 1000 K is a dopte d 
for O + CO surface reaction (Hasegawa & Herbst 1993). 



2.3. Chemical networks 

Three chemical networks are considered. The first one 
is the pure gas-phase network GAS. It consists of elec- 
trons, 12 atoms, 137 molecules, and 245 ions (395 species 
in total) involved in 3864 gas-phase reactions. The gas- 
grain network GAD is supplemented by 148 surface species 
and 729 gas-grain interaction processes (accretion, des- 
orption, and surface ion recombination). The gas-grain 
network DUST contains 192 additional surface reactions. 
These networks form the basis of our investigations. In 
the rest of the text we refer to them as "full" networks. 
Equations of chemical kinetics are solved with the stan- 
dard DVODE solver. Adopted solver parameters arc given 
in the Appendix. 



3. Reduction methods 

To reduce the number of species and/or reactions in a 
chemical network, one has to provide a mathematical 
method designed to estimate the significance of partic- 
ular species or reactions to the evolution of the species 
under consideration. In this paper we elaborate two such 
methods that arc outlined in this Section. 
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3.1. Species-based reduction 

This technique is derived from the method used by 
RRPHH and RBHPR. They have shown that it is pos- 
sible to reduce the number of species in a chemical net- 
work in order to compute abundances of selected species 
with a reasonable accuracy. The species-based reduction 
technique is based on the strategy to specify important 
species and then to select from the entire set only neces- 
sary species, that must be included in the chemical net- 
work in order to compute abundances of important species 
with a reasonable accuracy. 

Initially, the reduced species set consists of No impor- 
tant species only. For each species in the entire set we 
compute the value 



B l= £ 

3=1,1*0 



nj dfj 
9j dm 



where n, is the abundance of ith species, 
fj = hm,nin m - rij ^ kjini = Gj - L } 



(G) 



(7) 



is the net rate of the jth species abundance change which 
can be expressed as the difference between net gain Gj 
and loss Lj rates, and 



9] = max(Gj,Lj). 

With these definitions, the quantity 

rij dfj 
9] dm 



(8) 



(9) 



is the total rate of the abundance change of species j due 
to those reactions that involve species i only. The Bi is 
a measure of the sensitivity of the selected species abun- 
dances to the abundance of the ith species. As RRPHH 
noted, in this set of Bi values there is a clear boundary 
between "needed" and "not-so-needed" species. This is il- 
lustrated in Fig. |], where representative Bi values and a 
cut-off threshold t\ are shown. We select only those species 
from the entire set that have Bi above a boundary value, 
that is 



Bi > ei 



(10) 



They are added to the reduced set if they are not there 
already, and B calculation is repeated with the expanded 
species set. The process goes on until at some step no 
species are added to the reduced set. The entire network 
then consists of these species and of only those reactions 
that involve the selected species. Normally, the resultant 
reduced set is a combination of several intermediate sets 
constructed at different times to account for "early" chem- 
istry, "late" chemistry etc. But it turned out that the set 
reduced at the last time step (t = 10 7 years) works quite 
well at earlier times as well. 

It should be kept in mind that this method represents 
only the first part of the algorithm proposed by RRPHH. 




80 120 160 200 
Species Index 

Fig. 1. Representative B values for the HM-GAS-DIFF 
model computed at t = 10 7 years. Important species are 
Mg + and S + . The value of the cut-off parameter e\ is 
indicated. 
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-10 







1 00 200 
Reaction Index 



300 



Fig. 2. Representative w r values for the HM-GAS-DENS 
model computed at t = 10 7 years. An important species 
is CO. The value of the cut-off parameter €2 is indicated. 



3.2. Reaction-based reduction 

This new method allows to select only those species and 
reactions from the entire network that are necessary to 
compute abundances of important species with a reason- 
able accuracy. Unlike the species-based reduction tech- 
nique, where the number of reactions is decreased only 
as a by-product of the lower number of species, in this 
method the numbers of unnecessary species and ineffec- 
tive reactions are reduced simultaneously. In addition to 
the significantly shortened computational time, this could 
be useful if one aims at key destruction and formation 
pathways for important species, where a smaller number 
of chemical reactions makes the analysis easier. 

The basic idea of the reaction-based reduction is to 
search for the production and destruction reactions, most 
important for the evolution of an important species, and 
to determine their relative significance. This is done via 
the analysis how sensitive the net formation (or loss) rate 
of a given species is to the presence of particular reactions 
at a certain physical time. 
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At the beginning, the chemical model is computed 
with the full network in order to obtain abundances of 
all species from the network during the entire evolution- 
ary time span. 

Second, the important species are specified. Further, 
the algorithm estimates the weights of all species w s and 
reactions u> r in order to quantify their importance for the 
evolution of the species under investigation by the follow- 
ing iterative process. 

At the first iteration, weights of important species are 
set to 1 and weights of all other species are set to 0. At 
the ith iteration all possible formation and destruction 
pathways of the current species Sj are found and their 
significance values w\ are specified as 

<U) = 



Table 3. Reductions made for ionisation degree in a dif- 
fuse cloud. 



max < w r 



L (i), 



kjn ri (j)n r2 (j) 
^ hn ri {l)n r2 (l) 

/=l,AT r (i) 



w s {i) 



(11) 



Here kj is the rate of the jth reaction, n ri (j) and n r2 (j) 
are the abundances of the first r±(j) and second r2(i) reac- 
tants in the jth reaction, respectively, N T (i) is the number 
of reactions, in which Si is a reactant or a product, and 
w s (i) is the weight of the species Sj. 

Consequently, a new set of species, which are found at 
that iteration to be necessary for the evolution of the im- 
portant species and which were not considered at previous 
iterations, is formed. Their weights w s are estimated as 
the maximum possible values of the weights w r of the re- 
actions they are involved in. Then these species are added 
to the set of necessary species, and the last two steps are 
repeated. This process is performed for several stages dur- 
ing the computational time. 

The iterations are stopped when all species and all 
time moments are considered. Then one obtains a reduced 
chemical network from the full network by choosing only 
those reactions, that have weights u> r exceeding a pre- 
defined cut-off threshold €2 



logw r > e 2 . 



(12) 



The typical w r values and cut-off parameter 62 are 
shown in Fig. ||. The value of this cut-off is chosen to 
satisfy the requested accuracy requirement. In our calcu- 
lations the maximum allowed error is set to 30%. If there is 
a difference greater than 30% in abundances of the impor- 
tant species computed with the full and reduced networks 
for any particular time step, the cut-off is readjusted to a 
new, smaller value, and the last step is repeated. 

4. Ionisation degree 

We use the two reduction techniques described in the pre- 
vious section to build the reduced chemical networks that 
could be used for the computation of the ionisation de- 
gree in molecular clouds. All the models are first computed 



Model 


Important 


Reduced 


Speed 




species 


network 


gain 


Species-based reduction 


HM-GAS-DIFF 


S+ 


123 x 1583 


~ 25 


HM-GAS-DIFF 


S+, Mg+ 


126 x 1619 


~ 20 


LM-GAS-DIFF 




282 x 2539 


rv, 4 


LM-GAS-DIFF 


H+, HCO+ 


309 x 3177 


~ 3 


Reaction-based reduction 


HM-GAS-DIFF 


e~ 


58 x 111 


550 


LM-GAS-DIFF 


c~ 


73 x 169 


490 


HM-DUST-DIFF 


e~ 


210 x 683 


15 


LM-DUST-DIFF 


e~ 


357 x 1648 


2 



with the full chemical networks. Then we apply the reduc- 
tion techniques and compare the ionisation degree com- 
puted with the full and reduced networks. The results of 
the reduction together with relevant computational speed 
gains are summarised in Tables || ||. 

4.1. Diffuse cloud 

This case corresponds to the cloud at the initial stage 
of the protostellar collapse (nn = 10 3 cm -3 ). For pure 
gas-phase chemistry and high initial metal abundances 
(HM-GAS-DIFF), the dominant ions during most of the 
10 7 years are S + and Mg + . We found that it is suffi- 
cient to consider only these ions as important species in 
the species-based reduction method. Taking sulphur to be 
the only important species, we can reduce the number 
of species from 395 to 123 and achieve a computational 
speed gain of a factor of about 25 with uncertainty of 
~ 15% at t ~ 10 5 years (see Fig. ||, middle and bottom 
panels, squares). If we consider both ionised magnesium 
and ionised sulphur to be important species, the number of 
species in the reduced set is 126, and the accuracy is bet- 
ter than 5% during the entire evolution (see Fig. 0, middle 
and bottom panels, stars). The computational gain is al- 
most the same (see Table ||). Note that metals other than 
magnesium (i.e., sodium and iron) are almost equally im- 
portant for the ionisation degree. But as they enter the 
chemical network in a similar way, it is sufficient to des- 
ignate magnesium to be an important species in order to 
include in the reduced set all the species that are needed 
to compute abundances of Na and Fe as well. 

With the reaction-based reduction technique, we 
choose electrons to be the only important species. In 
this case the number of necessary species is 58, which 
is more than two times smaller compared to the species- 
based reduction. Additionally, the number of chemical re- 
actions is reduced from 3864 to 111 (see Fig. ||, top panel). 
The achieved computational time gain is of the order of 
500, while the resulting uncertainties are by definition 
smaller than 30% during all 10 7 years of the evolution 
(see Table 0). It should be recalled that the accuracy is 
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HM, full: 395 

HM, reduced (1): 123 

HM, reduced (2): 126 

LM, full: 395 

LM, reduced (1): 282 

LM, reduced (2): 309 



25 
20 
15 



JS 10 

CD 

rr 



I I 

□-oHM, reduced (1): 


123 


i 1 i 




-#■ *HM, reduced (2): 


126 






O-oLM, reduced (1): 


282 






x- xLM, reduced (2): 


309 










X 










if- 














I • • 


■'■x B 



Log 10 (Time, Years) 

Fig. 3. Ionisation degree (top and middle panels) and rel- 
ative errors (bottom panel) in a diffuse cloud, computed 
with the full and reduced networks for pure gas-phase 
chemistry (HM-GAS-DIFF and LM-GAS-D1FF models). 
Top panel reproduces results of the reaction-based reduc- 
tion. Indicated are the "exact" solutions for high (dotted 
line) and low (solid line) metallicities as well as the solu- 
tion for the reduced species sets (HM - dashed line, LM 
- long-dashed line). Middle and bottom panels illustrate 
the results of the species-based reduction. Shown are the 
"exact" solutions (HM - dashed line, LM - solid line) and 
the solutions for various reduced sets (see text). 

preset in the reaction-based reduction, so in all figures we 
only show relative errors for the species-based technique. 

The small number of necessary reactions in the re- 
duced network can be explained by the fact that in this 
case the chemistry is governed by a rather restricted set 
of ionisation-recombination processes. Such a set is easily 



extracted from the entire network by the reaction-based 
selection algorithm. The fundamental reactions during the 
entire evolution time are the ionisation of molecular hydro- 
gen and helium by cosmic rays and ionisation of carbon 
and sulphur by cosmic ray induced UV photons. These 
ions rapidly react with abundant molecules by efficient 
ion-molecule reactions, forming complex molecular ions, 
like HCO + . In turn, these molecular ions transfer charge to 
metal atoms. As metal ions are less chemically active, they 
become the most abundant charged species, regulating the 
fractional ionisation (Oppenheimer & Dalgarno 1974). 

The situation is more complicated from the chemical 
point of view when the initial metallicity is low. In the 
absence of abundant metals, the role of electron suppliers 
goes to H^, HCO + , and C + , involved in a rich chem- 
istry. Even though Hjj~ is the dominant ion, it alone does 
not determine the ionisation degree with enough accuracy 
which is illustrated in Fig. || (middle and bottom pan- 
els, circles). When this ion only is considered to be an 
important species in the species-based reduction, errors 
are negligible only at earlier times when the chemistry is 
relatively simple. However, after 10 4 years the errors grow 
significantly, exceeding a factor of 3 at the end of the com- 
putation. To account for later chemistry, one has to take 
into account another important ion, HCO + . When these 
two species are designated as important, errors do not ex- 
ceed 15% during the entire computational time (crosses 
in middle and bottom panels in Fig. ||). However, due to 
the complicated chemical connections of the two ions, the 
species set reduction is only modest in this case. Less than 
one hundred species can be excluded from the set with a 
reduction of the computational time by a factor of 3 only 
(see Table |). 

Not all the reactions involving these ions are impor- 
tant, though. With the reaction-based reduction method, 
we succeed in isolating a three time smaller set of 73 
species involved in 169 reactions which are needed to fol- 
low the evolution of the ionisation degree in this case. It 
decreases the time, needed for the computation, by a fac- 
tor of 500 (see Table |J). The slightly larger number of 
the necessary species and reactions in this case compared 
to the case of high metal initial abundances can be ex- 
plained by more diverse chemistry involving interactions 
of complex species. The accuracy of the computed ionisa- 
tion degree is again better than 30%. 

Inclusion of accretion and desorption processes as well 
as surface reactions into these models does not lead to 
any noticeable changes, primarily, because the density 
is too low to produce significant ice mantles. Results 
of the species-based method show that it is possible to 
use the same reduced set in HM-GAS-DIFF, HM-GAD- 
DIFF, and HM-DUST-DIFF models without significant 
loss of accuracy. The same is true for the low metallicity 
models. However, applying the reaction-based approach 
to the gas-grain chemical network including surface reac- 
tions (DUST-DIFF), we found only moderate reduction 
as nearly half of the species have to be retained in the 
reduced network. The reason is that the surface reactions 
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Table 4. Reductions made for the ionisation degree in a 



dense cloud with pure gas-phase chemistry. 


Model Important Reduced 
species network 


Speed 
gain 


Species-based reduction 


HM-GAS-DENS Na+ 72 x 568 
HM-GAS-DENS Na+,S+ 121 x 1560 
LM-GAS-DENS Na+ 71 x 567 
LM-GAS-DENS Na+,S+ 132 x 1717 


~ 100 
~ 30 
~ 100 
~ 25 


Reaction-based reduction 


HM-GAS-DENS e" 131 x 313 
LM-GAS-DENS e" 157 x 386 


130 

55 



have very high rates compared to the gas-phase reactions 
and due to the selectivity of the reaction-based method 
many of them are designated to remain in the reduced net- 
work. This leads to the larger number of chemical species 
and reactions in such a network compared to the pure 
gas-phase case and, consequently, to the smaller compu- 
tational time gains, which are about a factor of ten in this 
case (see Table |). 

4.2. Dense cloud 

4.2.1. Gas-phase chemistry 

When no gas-dust interaction is taken into account, the 
case of a dense cloud (nn = 10 7 cm' 3 ) at high metallicity 
is qualitatively similar to the low-density case. Metals are 
again the dominant ions. To compute their abundances 
(and hence the ionisation degree) with better than 50% 
mean accuracy one needs to include sodium only to the 
list of important species in the case of the species-based re- 
duction. The reduced set contains 72 species and provides 
two orders of magnitude gain in computational speed. To 
further improve accuracy we add sulphur to the list of 
important species. With the reduced set of 121 species 
we obtain errors less than 10% during the entire compu- 
tational time. The reaction-based technique gives similar 
results and produces a reduced network consisting of 131 
species and 313 reactions with the same computational 
gain. The difference between the high- and low-metallicity 
cases is smaller for a dense cloud than for a diffuse cloud. 
In a denser environment such chemically active ions as H;j" 
and HCO + are less abundant, and metal ions dominate 
the ionisation degree both in high- and low-metallicity 
cases. As a result of this similarity, in a dense cloud for 
the high and low metallicities the reduced species sets are 
nearly the same both for the species- and reaction-based 
approaches, with similar errors and computational gains 
(see Table @). 

4.3. Gas-grain chemistry 

The situation changes drastically, when gas-grain interac- 
tions are taken into account. In the high-metallicity case, 
at early stages of the evolution, the dominant ions are the 



same as in the diffuse medium, i.e., metals. However, af- 
ter 10 4 years their depletion becomes important, and the 
dominant ions are H+, HCO+, and N 2 H+. The straight- 
forward use of the species-based criterion under such con- 
ditions leads to the significant loss of accuracy at this time. 
The problem is that this criterion is based on the selec- 
tion of species whose abundances are needed to compute 
the abundance of an important species. But the depletion 
case is not that simple. Abundances of mantle species with 
high desorption energies are definitely not important and 
do not have to be computed with any accuracy. Based on 
the formal criterion (|l0|) , these species are to be excluded 
from the reduced set. However, this leads to the exclusion 
of the corresponding accretion reactions from the network, 
and counterparts of the excluded mantle species remain in 
the gas phase, significantly altering its chemistry. 

With the species-based reduction method we investi- 
gated the GAD-DENS model, i.e., the model in which ac- 
cretion and desorption processes are taken into account 
in addition to gas-phase chemistry. As we have no formal 
way to select species whose absence in the gas phase is 
important to compute the ionisation degree, we choose to 
check two artificial possibilities. First, we add surface met- 
als to the reduced network as they are primary contribu- 
tors to the electron abundance in the gas phase. Second, 
we allow all the selected neutral species to stick to dust 
grains, reasoning that if their gas phase abundances are 
deemed important then their depletion is important as 
well. Results are shown in Fig. || (middle and bottom pan- 
els). When only metal depletion (with a handful of other 
species selected automatically) is taken into account the 
errors are quite large, sometimes reaching 50%. To de- 
crease errors, depletion of all selected neutral species is 
needed. In the latter case, the approximate solution dif- 
fers from the "exact" one by as much as 20% but later on 
the difference is much smaller dropping to just a few per 
cent at t = 10 years (see Table |). 

The ambiguity of the previous selection procedure is 
not an issue in the reaction-based approach which is specif- 
ically designed to find necessary formation and destruction 
pathways rather than necessary species. As it was men- 
tioned above, this method of reduction treats gas-grain 
interactions and surface reactions equally with gas-phase 
reactions. Thus, one does not have to worry about which 
surface species to keep in a reduced set in order to esti- 
mate the abundances of important species. Therefore, this 
method works much better in the reduction of chemical 
networks involving complicated processes like gas-grains 
interactions and surface reactions. With the reaction- 
based technique, we reduced the number of species and 
reactions in the gas-grain chemical network involving sur- 
face reactions (DUST-DENS). We were able to remove 
about half of the species from the full gas-grain network 
and to reduce the number of reactions by a factor of 5 (see 
Fig. ^, top panel). The reduction speeds up the computa- 
tion by a factor of 8 only (see Table ||). At times, earlier 
than about 10 4 years, the ionisation degree is calculated 
with an accuracy better than 10%. At t > 10 4 years, the 
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Table 5. Reductions made for the ionisation degree in a 
dense cloud with gas-grain chemistry included. 



Model 


Important 


Reduced 


Speed 




species 


network 


gain 


Species-based reduction 


HM-GAD-DENS 




323 x 3201 


~ 6 




surface metals 






HM-GAD-DENS 


H 3 + , 


445 x 3445 


~ 2 




all surface 






LM-GAD-DENS 




325 x 2933 


r~j 4 




surface metals 






LM-GAD-DENS 


H 3 , 


444 x 3171 


~ 2.5 




all surface 






Reaction-based reduction 


HM-DUST-DENS 


e 


298 x 1060 


8 


LM-DUST-DENS 


e~ 


291 x 924 


7 



uncertainties grow but remain less than 30%. The rea- 
son is that it is more difficult to follow the evolution of 
the ionisation degree with the reduced network, since the 
importance of the complicated gas-grain interactions and 
surface reactions is raising up. In the low-metallicity case 
(LM-GAD-DENS), which is more complex from the chem- 
ical point of view, the reaction-based approach produces 
a reduced network that consists of 291 species and 924 re- 
actions and predicts an ionisation degree with a requested 
accuracy of 30% during the entire evolutionary time. It ac- 
celerates the computations by a factor of 7 (see Table . 

Comparison of top (DUST) and middle (GAD) panels 
in Fig. [| shows that there is a difference between the evo- 
lution of ionisation degree in the case of gas-grain chem- 
istry with and without surface reactions. When surface 
reactions are accounted for, fractional ionisation drops to 
10~ 9 in t ~ 10 4 years. In the case of gas-phase chemistry 
with accretion and desorption it reaches the same equilib- 
rium value later, after 10 5 years. Also, the maximum value 
of the ionisation degree is higher in the latter case. Thus, 
we conclude that surface reactions are important for the 
accurate estimation of the fractional ionisation in dense 
molecular clouds at late stages of their evolution. 

To conclude this section, we emphasise that the 
reaction-based technique is more suitable for the reduc- 
tion of gas-phase or gas-grain chemical networks if one 
intends to follow the evolution of the ionisation degree in 
diffuse or dense molecular clouds with a relatively small 
number of chemical species and reactions and a reasonable 
accuracy than it is the case for the species-based scheme. 

5. Carbon monoxide 

Essentially the same scheme of reduction, as in the previ- 
ous section, is utilised for CO as well. First, the chemical 
model is computed with one of the full chemical networks. 
Its results are used to select key species needed to predict 
the CO abundance. Then the validity of a reduced network 
is checked by comparing the abundances obtained with the 
full and reduced networks for the entire evolutionary time. 



In all cases, we assumed that the CO molecule is the only 
important species. Results of the reduction and relevant 
computational speed gains are presented in Tables ^ ||. 

5.1. Diffuse cloud 

In the case of pure gas-phase chemistry in a diffuse cloud 
(at the gas density jih = 10 3 cm~ 3 , GAS-DIFF) the main 
processes, relevant to the evolution of carbon monoxide, 
are its destruction by cosmic rays or by cosmic ray in- 
duced photons and its production through the neutral- 
neutral reactions of oxygen with carbon atoms and vari- 
ous light carbon-bearing molecules, dissociative recombi- 
nation of HCO + , and destruction of H^CO by CR-induced 
photons. Since metals as well as S, Si, P, and CI do not 
affect the evolution of CO abundances directly, it is nat- 
ural to expect that there should not be large differences 
between the high- and low- metallicity cases. 

In the absence of gas-grain interactions the evolution 
of CO is comparatively simple - during almost the entire 
evolutionary time its abundance is steadily approaching 
the equilibrium value of about 10 -4 which is reached at 
~ 10 6 years (see Fig. ||, top and middle panels). 

With the species-based reduction, of all 395 species we 
select only 36 species in the high-met allicity case and 27 
species in the case of low metal abundances (see Fig. |J 
middle panel) . Due to the small size of these reduced net- 
works, time needed to solve the relevant ODE system is 
negligible, less than a few seconds (see Table At early 
evolutionary phase (t ~ 10 3 — 10 4 years) the uncertainties 
of the CO abundance can reach 50% for high metallic- 
ity and 80% for low metallicity. Later, when the equilib- 
rium abundance is reached, error drops to about 20% (see 
Fig. ||, bottom panel). The reason for the high uncertainty 
at early times is that such rare species as CH3CO + and 
CH 2 CO, which are not recognised by the species-based al- 
gorithm to be necessary species, are in fact important at 
the early evolutionary phase for the accurate estimation 
of the CO abundance. 

Applying the reaction-based reduction, we are able to 
follow the evolution of the CO abundance for 10 7 years 
with an accuracy better than 30%. However, the price to 
pay for this is a significant number of species to be held 
in the corresponding reduced networks. One needs to keep 
158 and 95 species in the reduced networks for high and 
low metallicities, respectively (see Fig. [|, top panel). The 
relevant computational speed gain is a factor of 50 for the 
HM case and about 200 for the LM case (see Table |). It 
can be seen from Fig. || (top and middle panels), that the 
CO abundances in the case of high- and low-metallicities 
do behave similarly. Nevertheless, from the chemical point 
of view the case of low initial abundances of metals is 
somewhat simpler in comparison with the high metallicity 
case and implies smaller number of species and reactions 
to be kept in the reduced network. 

Inclusion of accretion-desorption processes and surface 
reactions does not affect the chemistry in a diffuse cloud 
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Table 6. Reductions made for CO in a diffuse cloud. 



Model 


Important 


Reduced 


Speed 




species 


network 


gain 


Spec 


les-based reduction 




HM-GAS-DIFF 


CO 


36 x 240 


> 100 


LM-GAS-DIFF 


CO 


28 x 182 


> 100 


HM-GAD-DIFF 


CO 


37 x 242 


> 100 


LM-GAD-DIFF 


CO 


27 x 173 


> 100 


Reaction-based reduction 


HM-GAS-DIFF 


CO 


158 x 588 


50 


LM-GAS-DIFF 


CO 


95 x 313 


240 


HM-DUST-DIFF 


CO 


298 x 1293 


8 


LM-DUST-DIFF 


CO 


306 x 1310 


5 



significantly. As we already mentioned in Sect. 4.1, the 
gas-grain interactions are not efficient at such low density 
(riH = 10 3 cm -3 ) since the typical rate of collisions be- 
tween dust grains and gas species is far too low. Using the 
species-based approach, we found that it is enough to re- 
tain only 37 species for the high metal initial abundances 
and 27 species in the LM case (see Table |). These arc 
almost the same numbers as in the case of pure gas-phase 
chemistry. In fact, one can use the same reduced sets of 
species and reactions in the GAS-DIFF, GAD-DIFF, and 
DUST-DIFF models (separately for each metallicity set). 
However, uncertainties in the computed CO abundances 
in all cases exceed 50% at t ~ 10 3 — 10 4 years. 

It is not that easy to remove these uncertainties. The 
application of the reaction-based method allows to cut 
the entire 543 species set by a factor of 2 in this case 
(see Table |). It is more than two times the number of 
species found to be necessary for the correct reproduction 
of CO abundances in the case of pure gas-phase chemistry. 
Consequently, it speeds up the calculations by a factor of 5 
only. The reason is, as noted in the previous section, that 
the reaction-based algorithm tends to keep some surface 
reactions in the reduced networks. Obviously, it leads to 
the larger number of species to be retained in the relevant 
reduced networks compared to the case of pure gas-phase 
chemistry. 

5.2. Dense cloud 

5.2.1. Gas-phase chemistry 

In a dense molecular cloud (nH = 10 7 cm -3 ) the main 
processes leading to the production and destruction of CO 
molecules in the case of pure gas-phase chemistry (GAS- 
DENS) are the same as in a diffuse cloud, but with one 
important change. At such high density, the role of chem- 
ical reactions involving abundant species like C, O, CH, 
and CH2, becomes even more important. Therefore, we 
may expect that in this case the chemistry relevant to 
the evolution of carbon monoxide is even simpler than in 
the diffuse cloud (GAS-DIFF). It implies a smaller set of 
necessary species. 



The results of reduction for this case are shown in 
Fig. |. The CO abundance increases rapidly at early evo- 
lutionary stages (t 55 10 3 years). Later, when the chemical 
equilibrium is reached, it remains constant around 10~ 4 . 

Using the species-based reduction technique, we suc- 
ceeded in the isolation of a small set of 27 species allowing 
us to follow the evolution of CO with a very good ac- 
curacy < 15% for all 10 7 years for both the high- and 
low-metallicity case (see Fig. ^, middle and bottom pan- 
els). Again, the computational time is negligibly small, 
less than a second (see Table [7]) . 

But even this 27 species set contains excessive informa- 
tion. Applying the reaction-based reduction approach, we 
revealed that in reality there are only 8 species involved 
in 9 reactions which are necessary for the accurate esti- 
mation of CO abundances during the entire evolutionary 
time (see Fig. ||, top panel). The relevant reduced net- 
work is given in Table ||. The uncertainties of computed 
CO abundances do not exceed 15% for all 10 7 years of 
the evolution. This remarkable result, showing the power 
of the reaction-based technique, explains why there is no 
difference between the cases of high and low metals in 
a dense cloud: metals are completely excluded from this 
network. 

In the medium with the atomic initial composition (ex- 
cept for molecular hydrogen) reactions start from the dis- 
sociation of H2 molecules by cosmic ray particles and from 
their radiative association with atomic carbon. The lat- 
ter reaction leads to CH2 formation. This molecule reacts 
with atomic H and C leading to CH and C2H. These three 
molecules (CH, CH2, and C2H) associate with atomic oxy- 
gen finally producing the CO molecule. These processes 
are balanced by destruction of CO and CH by CR-induced 
photons. In this way, the key reaction pathways can be 
isolated for which high-quality experimental reaction rate 
measurements have to be provided for an accurate mod- 
elling. It may seem to be somewhat surprising that the re- 
duced network does not contain HCO + , at variance with 
the diffuse cloud case. Note, however, that the maximum 
value of the HCO + abundance reached in the high-density 
model (~ 10~ 10 ) is two orders of magnitude less than in 
the diffuse cloud model. Given the correspondingly low 
electron abundance, this prevents effective CO synthesis in 
dissociative recombination of HCO + . Proton transfer re- 
actions of CO with H^ and H^" that would transform CO 
to HCO + are also nearly quenched due to low abundances 
of these ions. Thus, while CO is definitely an important 
species for HCO + in a dense medium, the opposite is not 
true. 

Going somewhat out of the scope of the present paper, 
we mention briefly that this network proved to be valid un- 
der a much wider range of conditions that are described 
here. Calculations that will be presented in a subsequent 
paper show that it can be used for estimating the CO 
abundance not only in dark dense clouds, but in a translu- 
cent (Ay > 1 mag) medium as well. It gives less than 50% 
error at 10 4 < n H < 10 10 cm~ 3 and T < 250 K and can 
be used in a denser medium as well (riH < 10 12 cm~ 3 ) 
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Table 7. Reductions made for CO in a dense cloud with 
pure gas-phase chemistry. 



Model 



Important Reduced Speed 
species network gain 



Species-based reduction 



HM-GAS-DENS CO 
LM-GAS-DENS CO 



27 x 133 > 100 
27 x 152 > 100 



Reaction-based reduction 



HM-GAS-DENS CO 
LM-GAS-DENS CO 



8x9 
8x9 



10 4 
10 4 



Table 8. The reduced gas-phase network for CO chem- 
istry in a dense cloud obtained with the reaction-based 
method. 



Reactions 



1 H + CH 2 

2 H 2 + C 

3 C + CH 2 

4 O + CH 

5 + CH 2 

6 + C 2 H 

7 H 2 + CRP a 

8 CH + CRPHOT b 

9 CO + CRPHOT 



CH + H 2 

CH 2 + PHOTON 

C 2 H + H 

CO + H 

CO + H + H 

CO + CH 

H + H 

C + H 

c + o 



a Cosmic ray particle 
b CRP-induced photon 



provided T ~ 200 — 250 K and Ay > 5 mag, conditions 
being typical of hot cores or protoplanetary discs. 

The accuracy of this network is also insensitive to the 
variation of the adopted C/O ratio (0.43) within a factor of 
a few. The value of this ratio would be critical, though, in 
the attempt to construct a reduced network for carbon or 
oxygen separately. If carbon is less abundant than oxygen 
(as in the adopted abundance sets), all C atoms are locked 
in CO molecules, and the reduction is rather effective. In 
fact, it is possible to construct a single network of a few 
tens of species that reproduces abundances of C, O, and 
CO simultaneously with less than 100% uncertainty. If the 
C/O ratio is greater than 1, extra carbon is available for 
producing long carbon chains, and the effective reduction 
(with C as an important species) is not possible. 

5.2.2. Gas-grain chemistry 

The situation is no longer trivial when gas-grain interac- 
tions and surface reactions are taken into account (GAD- 
DENS and DUST-DENS). The influence of these pro- 
cesses on the evolution of abundances is twofold. First, 
the freeze-out of the gas-phase CO leads to the activation 
of certain chemical reaction chains that are suppressed in 
the presence of CO. Second, the presence of CO in icy 
mantles can drive some surface reactions that may even- 
tually result in the formation of complex organic species. 



In the absence of surface reactions the evolution of car- 
bon monoxide is simple (see Fig. ^, middle panel). First, 
its abundance grows steadily, at t Ss 5 TO 3 years reaching a 
value of ~ 5T0~ 5 . Then the gas-grain interaction becomes 
important, and the CO molecules heavily deplete from the 
gas phase, forming significant mantles. The chemical equi- 
librium is reached at 10 5 years with the gas-phase CO 
abundance of a few times 10 -8 . 

The species-based reduction technique isolates two sets 
to be used in the high- and low-metallicity cases of the 
GAD-DENS model. But in fact these sets are very simi- 
lar to each other, and one can use the same set contain- 
ing about 30 species with high- and low- initial mctal- 
licities both in pure gas-phase models and with accre- 
tion/desorption taken into account. The resulting uncer- 
tainties are smaller than 10% for both metallicity cases 
(see Fig. [j| bottom panel). Again, the computational time 
which is necessary to compute the evolution of CO is neg- 
ligibly small. 

The evolution of carbon monoxide in the situation, 
when surface reactions are considered, is shown in Fig. 0, 
top panel. The gas-phase CO abundance increases at early 
times and reaches the value of ~ 10~ 4 at t ~ 10 3 years. 
Then, the depletion of CO becomes efficient. At t ~ 
10 4 — 10 6 years the gas-phase CO abundance is deter- 
mined by the balance between accretion and desorption 
processes. Its value is almost the same as the equilibrium 
CO abundance in the GAD-DENS model, i.e., a few times 
10 -8 . After 10 6 years surface re-processing of CO takes 
over. It is converted to other species, mainly, CO2. At 
this stage the accretion of CO is no more balanced by its 
desorption, and after some time the new, much lower equi- 
librium value ~ 10~ n is established. Note that this be- 
haviour is at least partially caused by our simplified treat- 
ment of th e gas- dust interactions. We adopt the Hasegawa 
& Herbst ( 1993 ) value for CO binding energy which is es- 
timated for a Si02 surface. In reality, after the CO mantle 
is formed, the binding energy of CO and other molecules 
may become smaller, so that desorption is still competitive 
with surface re-processing. 

With the reaction-based reduction, we are able to se- 
lect sets of 388 species and 1563 reactions for the high- 
metallicity case and 270 species and 798 reactions for the 
low-metallicity case. The relevant computational gains are 
only about 3 and 16, respectively (see Table ||). The error 
of the predicted CO abundances is less than 10% during 
the entire evolutionary time. The number of species which 
are needed to estimate the CO abundance accurately in 
the high-metallicity case is again ~ 50% larger compared 
to the low-metallicity case. Reduced networks with smaller 
number of species are not usable in the calculations in this 
case. The reason is that the relevant ODE systems are so 
stiff for the DVODE solver that the computational times 
are much larger compared to the case of the full chemical 
network. 

Concluding this section, we emphasise that if surface 
reactions are considered in a chemical network, then its 
reduction for carbon monoxide is of no practical use as 
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Table 9. Reductions made for CO in a dense cloud with 
gas-grain chemistry included. 



Model Important Reduced 
species network 


Speed 
gain 


Species-based reduction 


HM-GAD-DENS CO 31 x 148 
LM-GAD-DENS CO 29 x 133 


> 100 

> 100 


Reaction-based reduction 


HM-DUST-DENS CO 388 x 1563 
LM-DUST-DENS CO 270 x 798 


3 
16 



the relevant computational speed gains are too small. 
However, if accretion and desorption of species are taken 
into account as the only gas-grain processes, then it is pos- 
sible to find a few tens in half a thousand species which 
must be kept in a reduced network in order to follow the 
evolution of CO without significant (<: 50 — 100%) errors. 

6. Discussion 

The purpose of the present paper is to check whether a 
typical astrochemical database contains redundant infor- 
mation that can be ignored in certain calculations regard- 
ing interstellar chemistry. The primary incentive for our 
work is the long time usually needed to complete a more 
or less detailed model of interstellar chemistry in dynam- 
ically evolving media. In many cases the evolution of only 
a few or even one species are of the main interest for an 
astronomer. The notable example is the degree of ionisa- 
tion in dense interstellar clouds which is regulated by the 
abundances of a few dominant ions. It has been shown 
by many authors that in this case the relevant chem- 
istry is governed by a rather restricted set o f chemical 
processes (e. g. Op penheimer & Dalgarno 1974 ; Ciolek & 
Mouschovias 1995 ). 

Of course, one needs a mathematically rigorous way to 
"extract" only those species and reactions from the entire 
network that are needed to reproduce the evolution of a 
species under investigation within a reasonable numerical 
accuracy (say, < 100%). Such a reduction would decrease 
the number of species in the network and, hence, the rel- 
evant computational time. Using the DVODE solver, we 
found that this time is proportional to roughly the third 
power of the number of the involved chemical species. 
Thus, by decreasing the number of species in the model 
by just a factor of two, one can obtain an order of mag- 
nitude computational speed gain! Besides, the reduction 
would help to understand the role of a particular chemical 
compound or reaction for the evolution of an important 
species. However, one must bear in mind that the cost of 
the reduction is that the relevant reduced network has a 
certain range of applicability, i.e. it gives reliable estimates 
on the abundances of an important species for a specific 
range of the gas temperature, density, etc. 

The first mathematically proven reduction of astro- 
chemical networks is made by RRPHH. Their intention 
was to isolate a set of species and reactions which are 



needed to follow the evolution of the CO abundance within 
the 30% accuracy. They considered physical conditions, 
typical of translucent regions, i.e., nn = 500 — 10 5 cm~ 3 , 
T = 10 K, Ay = 0.5 — 5 mag, with "low metal" ini- 
tial abundances. They discriminated 69 necessary species 
participating in 241 reactions from the Ohio State New 
Standard Model, including reactions of dissociative recom- 
bination on dust grains. They found that it is possible to 
decrease further the number of necessary species and re- 
actions to 33 and 116, respectively, by the artificial mod- 
ification of the network. It consisted in exclusion of all 
reactions involving carbon chains longer than C2H2" and 
in introduction instead of few synthetic reactions to cor- 
rectly predict the abundances of shorter carbon chains. 

The same method was applied by RBHPR in order to 
identify reduced networks that predict the degree of ioni- 
sation in interstellar clouds with less than 30% uncertainty 
under a wide range of physical conditions. They reduced 
two networks, namely, a "full" network of 218 species and 
2747 reactions and a "large" network of 79 species and 
1026 reactions. In the "full" network metals beside sodium 
were excluded, whereas the "large" network in addition 
excluded nitrogen and most species having more than two 
nuclei heavier than hydrogen. High and low metal ini- 
tial abundances were used. They considered the following 
physical conditions: T — 20 K, 7Jh = 500 — 10 5 cm -3 , 
A\ = 1 — 10 mag. The small combined reduced network 
having 28 species and 71 reactions was found for the case 
of low metallicity. In the high metallicity case similar net- 
work of 29 species and 65 reactions was isolated with the 
artificially simplified sulphur chemistry. 

We analysed the possibility of reducing the number of 
species and reactions in the UMIST 95 database by two 
different methods that we refer to as species-based and 
reaction-based techniques. Our primary goal is to com- 
pare the efficiency of these techniques under physical con- 
ditions typical of dark regions. We study the applicability 
of both methods in the case, when the complicated gas- 
grain interactions and the surface reactions are taken into 
account. Having in mind possible applications, we tried 
to make these methods as automatic as possible. Unlike 
RRPHH and RBHPR, we did not apply artificial modifi- 
cations to the reduced networks. While such modifications 
can be useful in certain cases, usually it is not a trivial task 
and can be fraught with difficulties when a wide range of 
physical conditions is considered. Thus, we left the reduc- 
tion entirely to the reduction techniques. 

Both techniques proved to be efficient enough to re- 
duce significantly gas-phase chemical networks needed to 
compute the ionisation degree or CO abundance in diffuse 
and dense molecular clouds with a reasonable accuracy of 
< 30%. If the accretion-desorption processes and surface 
reactions are accounted for, the reduction becomes less 
effective since the relevant chemistry is too complicated. 
Even though the number of chemical reactions in the re- 
duced networks could be decreased by a few times, making 
their analysis easier, the computational speed gains are 
too small in these models to make this reduction of par- 
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ticular interest for chemo-dynamical modelling. Though, 
at least for the ionisation degree, the situation may be al- 
leviated by the fact that at high densities when depletion 
is effective, grain charging processes are more important 
for the fractional ionisation than the chemistry. 

We found that the reaction-based technique is more 
effective compared to the species-based technique (see 
Tables ||-|5|). The number of reactions is at least three 
times smaller for the reaction-based reduction than for 
the species-based reduction which is natural to expect as 
reaction rates are not discriminated in the latter method. 
But the number of species is also smaller for the reaction- 
based technique in most models. When the gas-dust inter- 
action is taken into account, the reaction-based technique 
is even more effective than the species-based technique, 
as it provides the simultaneous reduction of the number 
of species and reactions and proper treatment of the gas- 
grain processes. Unfortunately, it is impossible to compare 
our results for the ionisation degree directly with the re- 
sults of RBHPR. One reason is that they considered mod- 
erate values of visual extinction Ay = 1 — 10 mag, typical 
of translucent clouds, whereas we dealt with obscured re- 
gions, where Ay > 10 mag. Another issue is that RBHPR 
started the reduction from a somewhat smaller "full" net- 
work of 218 species and 2747 reactions, while we consid- 
ered the network of 395 species involved in 3864 reactions. 
Nonetheless, the resulting reduced networks have compa- 
rable sizes, namely, 29 species and 65 reactions (RBHPR, 
see Table 4) and 58 species and 111 reactions (this work, 
Table 3) in the case of high metals. In the low metallic- 
ity case the reduced network by RBHPR has 28 species 
involved in 71 reactions, compared to our network of 73 
species and 169 reactions. 

In the case of reduction made for CO, the situation 
is a bit more diverse. In a diffuse cloud with pure gas- 
phase chemistry, the species-based reduction technique 
overcomes the reaction-based method, but at the expense 
of higher uncertainties (see Table || and Fig. ||). In a dense 
cloud with pure gas-phase chemistry, the reaction-based 
technique allows us to extract the reduced network which 
is very small (Table ||) but still functional under a very 
wide range of physical conditions. This network is smaller 
than the one found by RRPHH, though it has been ex- 
tracted from the larger full network. Also, it has a wider 
range of applicability. 

The next question to ask is how to apply these methods 
in dynamical calculations. One can envisage two different 
strategies to pursue this goal. One way is to reduce the 
number of species and reactions at the outset and then to 
use the reduced network during the entire computational 
time. Another strategy would be to perform the network 
analysis in the "real-time" mode, during the dynamical 
evolution. 

The first strategy is obviously more cost-effective as 
it does not add any overhead to the computational time. 
The reduced network must be prepared in advance and 
checked against the range of physical conditions that are 
expected in the dynamical model. Of course, this check 



is only possible under simplified conditions (otherwise it 
makes no sense). For example, if the range of densities 
(temperatures) is expected in the model, one may check 
that the reduced network is valid at its upper and lower 
boundaries, hoping that it is valid as well for any den- 
sities (temperatures) in between. Of our two techniques, 
the reaction-based approach works in this mode, as to con- 
struct the reduced network the abundances during the en- 
tire evolutionary time are analysed. The resultant network 
accounts for both the early and later chemistry. 

The second method seems to be more robust. After 
the ability of the reduction technique to remove the ex- 
cessive information from the chemical network is checked, 
one can perform the reduction, for example, at each dy- 
namical time step or even less often. The full chemical 
network must be taken into account every time, with cur- 
rent abundances, temperature, density, etc. One may want 
to identify not only species that can be ignored but species 
with quasi-stationary abundances as well. Our species- 
based reduction technique is well suited to be used at a 
given t with the current values of parameters because, as 
we mentioned above, the sharp jump in Bi values allow 
to distinguish necessary and unnecessary species clearly. 
Usually, there is no clear boundary between needed and 
unneeded reactions (see Fig. 2), so one has to readjust 
cut-off threshold several times before optimal reduction 
is found. Thus, in some cases, the reaction-based method 
is more time-consuming than the species-based approach 
which is critical if the method is to be invoked from the 
dynamical code more than once. 

Neither method guarantees that the same accuracy, as 
found in test cases, is preserved during the real computa- 
tions. With the first method one cannot be sure that con- 
ditions, under which the reduced network is not valid, are 
never met in the dynamical model. The second method 
properly takes into account the current conditions but 
is more sensitive to the history of abundance changes. 
Consider the ionisation degree in the dense cloud when 
freeze-out is taken into account (high-metallicity case). 
One of the dominant ions at later times is N2H" 1 ". As a 
necessary species, it is included in the network at the time 
when metal depletion becomes important (late stages). 
However, the abundance of this ion cannot be followed 
with any accuracy as it is absent in early-time reduced 
networks, and its entire previous evolution is lost. So, both 
strategies must be applied with care, taking into account 
possible uncertainties. 

We note in passing that the decreased accuracy of the 
reduced networks may not be such an issue. Many reac- 
tions in the UMIST 95 have poorly measured or estimated 
rates, with typical uncertainties ^ 25 — 100% or more 
(Millar et al. 1997). This introduces unknown global er- 
rors in the modelled species abundances. On the other 
hand, modern observations usually provide column den- 
sities of abundant species with uncertainties of few tens 
of percents. Therefore, at nearly all astrophysically im- 
portant situations accuracy of ~ 50% seems to be quite 
appropriate. 
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7. Conclusions 

We developed a robust method to reduce simultaneously 
the number of species and reactions in chemical networks. 
The straightforward application for such reduced chemical 
networks would be the modelling of the evolution of mag- 
netised protostcllar clouds or protoplanetary discs, when 
it is necessary to compute the fractional ionisation self- 
consistently with dynamical processes. 

Investigating the applicability of this reduction ap- 
proach, we considered the conditions of diffuse and dense 
molecular clouds and utilised the UMIST95 database of 
chemical reaction rates. For the sake of comparison, we 
performed the same reduction with the objective reduc- 
tion technique. The new reaction-based way of reduction 
proved to be more efficient and accurate but more time- 
consuming as well. 

It is found that the number of species and reactions 
that are needed to follow the evolution of the CO abun- 
dance or the ionisation degree can be greatly decreased 
in the case of pure gas-phase chemistry, implying com- 
putational speed gains of more than a few hundred. For 
instance, to follow accurately (with < 15% uncertainty) 
the evolution of the CO abundance in a dense cloud, it is 
enough to retain only 8 species involved in 9 reactions in 
the chemical network. However, as soon as the gas-grain 
interactions and surface chemistry arc taken into account, 
reduction is modest and corresponding acceleration fac- 
tors are around 10. 
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Appendix A: Adopted input parameters 

We use the standard DVODE package designed to solve 
systems of ordinary differential equations. The relevant 
input parameters adopted in our computations are the 
following: 

— absolute error e abs = 2.20 • 10~ 16 ; 

— relative error e rcl = 10~ 8 ; 

— Jacobian computed by the solver (internally). 

We found that the relative error e re i = 10 -8 is enough to 
give reliable estimations of the abundances of the species. 

We use the g77 compiler and SuSe Linux 7.3 installed 
on the Pentium IV 2 GHz, 1 GB RAM PC as well as 
Compaq Visual Fortran 6.5 compiler on Pentium III, 
850 MHz, 256 MB RAM PC. The typical time necessary 
to run a model is only a few minutes in the case of gas- 
phase chemistry and can reach half an hour if the gas-grain 
interactions are taken into account. 
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Fig. 4. Ionisation degree (top and middle panels) and rel- 
ative errors (bottom panel) in a dense cloud, computed 
with the full and reduced networks for gas-phase chem- 
istry including accretion and desorption processes and 
dust surface reactions. Top panel reproduces results of 
the reaction-based reduction for the DUST-DENS model. 
Indicated are the "exact" solutions for high (dotted line) 
and low (solid line) metallicities as well as the solution for 
the reduced species sets (HM - dashed line, LM - long- 
dashed line). Middle and bottom panels illustrate results 
of the species-based reduction for the GAD-DENS model. 
Shown are the "exact" solutions (HM - dashed line, LM - 
solid line) and the solutions for various reduced sets (see 
text). 



Fig. 5. The CO abundances (top and middle panels) and 
relative errors (bottom panel) in a diffuse cloud, computed 
with the full and reduced species set for pure gas-phase 
chemistry (GAS-DIFF model). Top panel reproduces re- 
sults of the reaction-based reduction. Indicated are the 
"exact" solutions (HM - dotted line, LM - solid line) 
and the solution for the case when CO is considered to 
be the only important species (HM - dashed line, LM - 
long-dashed line). The middle and bottom panels illustrate 
results of the species-based reduction. Shown are the "ex- 
act" solutions (HM - dashed line, LM - solid line) and 
the solution for the case when CO is an important species 
(HM - open squares, LM - open circles). 
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Fig. 6. Same as in Fig. || but for a dense cloud (GAS- 
DENS). 



Fig. 7. Same as in Fig. | but for the DUST-DENS (top 
panel) and GAD-DENS (middle and bottom panels) mod- 
els. 



